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Abstract. We study the accretion efficiency of neutron stars in long- 
period binaries {P ^ 20^^) wliich accrete from a giant companion. Using 
a-disc models and taking into account the effect of irradiation of the 
accretion disc by the central accretion light source we derive explicit ex- 
pressions for the duty cycle d and the accretion efficiency r] in terms of 
the parameters of the binary system and the disc instability limit cycle. 
We show that the absence of millisecond pulsars in wide binaries with 
circular orbits and periods P ^ 200*^ can be understood as a consequence 
of the disc instability if the duration of the quiescent phase between two 
subsequent outbursts is at least a few decades. 



1. Introduction 

Ever since the detection of the first millisecond pulsar (Backer et al. 1982) it 
has been clear that a ms-pulsar is a neutron star (NS) which has been spun up 
by accretion in a close binary system (for a review see Bhattacharya & van den 
Heuvel 1991). Because the amount of mass needed to spin up a NS to ms periods 
is rather small, i.e. O.O5M0 ^ AMns ^ O.IMq (Burderi et al. 1999) and the 
mass available from the donor is typically many times larger, i.e. typically of 
order IMq, it appeared as if there was no problem spinning up a NS to such 
short periods. Yet in recent years, evidence has accumulated which shows that 
NS are inefficient accretors. Thorsett & Chakrabarty (1999) noted that in all 
cases where the mass of a NS could be measured with precision it turned out 
to be in the very narrow interval I.SMq ^ Mns ^ 1.4M0. Even if one makes 
the extreme assumption that all NSs are born with exactly the same mass, this 
would imply that none of these NSs could have accreted more than ~ O.IM0, 
despite the fact that some of these NSs have been spun up by accretion and the 
mass available from the donor star was probably much larger than ~ O.IMq. 
In another context. King & Ritter (1999) have presented further evidence that 
NSs can be very inefficient accretors. The main reason for this low accretion 
efficiency is the frequent disparity between the mass transfer rate provided by 
the donor star and the maximum rate at which the NS can accrete. The latter is 
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approximately the Eddington rate Me ~ 1.5 x lO~^M0yr~^, whereas, depending 
on circumstances, the former can be larger by many orders of magnitude. 

One of the obvious places in which NSs can be spun up to ms-pulsars are 
wide binaries in which a NS accretes from a giant. Tauris &: Savonije (1999) 
have recently devoted a detailed study to this formation channel. One of their 
predictions is that, in principle, ms-pulsars should form in binaries with final 
orbital periods of up to P ~ lOOC^, even when taking into account that the 
mass transfer rate in long-period systems (P ^ 20^^) exceeds the Eddington 
accretion rate. In addition, their calculations predict a correlation between the 
final orbital period and the mass of the ms-pulsar. Unfortunately, this cannot 
be tested at present. But their first prediction can. And here observations show, 
contrary to their prediction, that there are no ms-pulsars in wide binaries with 
a circular orbit and a period P ^ 200^^ (see e.g. table 1 in Taam, King Sz Ritter 
2000). In this context, Li & Wang (1998) have already noted that the accretion 
efficiency of the NS in such a binary is far lower than has been assumed by 
Tauris &: Savonije (1999) if one takes into account the fact that in all these wide 
binaries the accretion disc around the NS is thermally unstable and accretion 
is transient (King et al. 1997). This is because a) during an outburst the mass 
flow rate through the disc is much larger than the mass transfer rate from the 
donor, which, in turn, is already supcr-Eddington if P 20^^, and b) because a 
NS cannot accrete in excess of the Eddington rate. In fact, Li &; Wang (1998) 
found that if the duty cycle of the disc instability limit cycle is d ^ 10~^ then 
in systems with a final orbital period Pf ^ IOC* the total mass accreted by the 
NS is AMns ^ few 1O~^M0, too little to spin up a NS to ms periods. 

In this paper we elaborate on the idea of Li & Wang (1998) and derive 
explicit expressions for the duty cycle and the accretion efficiency of the NS in 
terms of the parameters of the binary system and the outburst cycle. Our calcu- 
lations show that if the duration of the quiescent phase between two subsequent 
outbursts exceeds a few decades, the accretion efficiency becomes so small that 
in binaries with a final orbital period Pf ^ 200'*, the NS cannot accrete enough 
mass to become a ms pulsar. 



2. Basic model ingredients 

The binary systems for which we wish to decribe the accretion processes are 
wide, long-period binaries with orbital periods of typically P ^ 20^* and in 
which a neutron star primary of mass Mi accretes from a giant secondary, of 
mass M2. The basic model ingredients we use can be summarized as follows: 

2.1. Mass transfer 

We assume that mass transfer from the secondary is thermally (and adiabati- 
cally) stable. In this case, the mass transfer rate can be computed from a simple 
analytical approximation (e.g. Ritter 1999). Following King et al. (1997) and 
using the same model parameters we can write for the mass transfer rate 

-M, = 7.3 X 10-^0Meyr-i(Ce - CR)-'ml-'''V-'''\ (1) 
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where 1712 = M2/M.Q,p = P{d) is the orbital period in days, the thermal 
equilibrium mass-radius exponent of the donor star, and the mass-radius 
exponent of the donor's critical Roche radius. For (thermal) stability we require 
Ce - Cil > (e-g- Ritter 1988). 

2.2. Disc irradiation 

Accretion on to a compact star can result in significant external irradiation of 
the accretion disc by the central accretion source (e.g. van Paradijs 1996). Here 
we follow King & Ritter (1998) and assume that the effective temperature of 
the disc is kept above the hydrogen ionization temperature Th 6500K by 
irradiation from the central accretion source out to a radius 

i?h = (SlMe)^/^ (2) 

where Mg is the central accretion rate, and, for a neutron star accretor, Bi = 
3.9 X 10^ cgs (King &; Ritter 1998). Furthermore, we shall assume that the 
neutron star cannot accrete in excess of the Eddington accretion rate Me — 
1.5 X lO~^M0yr~^, i.e. that Mc ^ Me- This means that the maximum radius 
out to which the disc can be kept hot by irradiation, i.e. at Tes > Th, is 

i?crit =^7.1R0(&imE)'/'. (3) 

Here ttje = Me/W^^Mqyt^^ and 61 = Bi/3.9 x 10^ cgs. At the same time our 
assumption that Mc ^ Me means that whenever the mass flow rate through the 
disc is — Mfj > Me, the excess matter must be ejected from the disc. 

2.3. The a disc model 

For describing the properties of the accretion disc around the neutron star we 
use the Shakura Sz Sunyaev (1973) a-disc model. Accordingly, if the vertical 
temperature stratification of the disc is dominated by viscous heating, the disc 
viscosity is (e.g. Frank, King & Rainc 1992) 

i^visc = 2.7 X 10i5cm2s-^a4/5(-md)'/'Vr'^'^'/', (4) 

where — Md is the mass flow rate through the disc, R the disc radius, rhd = 
Md/10~^MQyr~^, and r = R/Rq. a is the viscosity parameter, here assumed 
to be constant with radius. 

(4) holds only to the extent that the vertical temperature stratification in 
the disc is dominated by viscous heating, hence its dependence on Mj. If, on 
the other hand, the disc is strongly heated by external irradiation, the vertical 
temperature profile becomes nearly isothermal at the irradiation temperature 
Tirr- Taking this into account in the a-viscosity ansatz we find (see also King 
1998) 



a^(BiMc)^/^(GMi)-V2i? 



= 3.8 X lOi^cm^s-^a (^)'^ {hm,y/^m^''\. 



Here TZ is the gas constant and ji the mean molecular weight. 
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The maximum suface density which a disc can reach in the cool state, 
i.e. before it becomes thermally unstable and must turn into the hot state, 
is conventionally denoted by Sb- Taking from computations of Ludwig, 
Meyer-Hofmeister & Ritter (1995), the mass accumulated before an outburst 
sets in can be written as 

M<i ~ 7 X lO-'^Me a-^/Vr'^-'V^-i'^J^. (6) 

Here = R,i/RQ, where is the outer radius of those parts of the disc which 
are involved in the outburst, ct^ the a-paramctcr in the cool state, and < 1 is 
the filling factor defined by this equation. takes into acount that a the onset 
of an outburst S = at only one point in the disc and that S < Sb elsewhere. 

2.4. Assumptions about the disc instability limit cycle 

We shall assume that the disc instability limit cycle in these wide, long-period 
binary systems with a neutron star accretor can be described as follows: 

a) During an outburst the disc is kept in the hot state by irradiation out to the 

radius Rh{Mc) given by (2). Because Mc ^ M^,R\,{Mc) ^ Rh{M^)- If 
the physical outer radius of the disc -Raise > Rh{^c)^ only the inner part 
with radius R^ = R^^Mc) is hot, whereas the outer part with i?h(-^c) < 
R < -Rdisc remains in the cold state (e.g. King & Ritter 1998). Hereafter 
we shall call R^ the radius of the "active" disc. If, on the other hand, 
Rh{Mc) > i?disc)-Rd = -Rdisc- If the mass flow rate through the outer 
parts of the active disc is —M^iRd) > Me, then as long as — Md(-Rd) > 
Me, Mc = Me, and Rh{Mc) = Rh{ME) = const. 

b) During an outburst the active part of the disc (out to R^) is essentially 

emptied, i.e. the mass remaining in the active part of the disc immediately 
after an outburst is insignificant compared to the mass M(j given by (6) 
immediately before the onset of an outburst. 

c) During quiescence the mass of the active part of disc which was accreted/ 

ejected during the previous outburst is replenished at the mass transfer 
rate — M2. This is tantamount to assuming that if an inactive outer part 
of the disc exists, i.e if i?h(-M^c) < -Rdisc, this outer part is stationary, i.e. 
that Md(-R) = M2 for all Rh{Mc) <R< -Rdisc- 

d) We assume that the hot part of the disc (during an outburst) is characterized 

by a viscosity parameter O'l^. whereas in the cool state (i.e. in the active 
part during quiescence and in the inactive part) the viscosity parameter is 
Q!c < ah- 

3. Disc accretion in long— period neutron star binaries 

Here we are mainly interested in systems which end their nuclear time scale- 
driven mass transfer at very long orbital periods, i.e. P ^ 200^^. With the 
analytical solution (Ritter 1999) one can show that the corresponding initial 
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periods at which stable mass transfer must have started are typically Pi ^ 20^'. 
(the precise value depending on whether mass transfer is conservative or not, cf. 
Taam et al. 2000). It has been shown by King et al. (1997) and Li & Wang 
(1999) that at such long periods the mass transfer rate provided by nuclear 
evolution of the donor is always too small for stable disc accretion, i.e. the 
systems in question are transient, going through a disc instability limit cycle. 

In addition, the orbital period at which -Rdisc = -Rcrit (with -Rdisc > ^crit for 
longer periods) is 



crit 



12d 



-^disc/ Rl,R 
0.7 



-3/2 



m 



-0•67V0■l^^(6lmE)'/^ 



(7) 



where Ri^r is the critical Roche radius of the primary. Thus, for practically all 
systems in question we also have i?disc > Rh{Mc)- 

As we shall argue below, during most of an outburst the mass flow rate 
through the outer parts of the active disc is highly super-Eddington in the sense 
that — Md(^d) ^ M-^. Therefore, at least initially, in these systems = 
Ry,{Me) = const. 

Now, if the viscosity at i? = Rd docs not depend on Md, and i?d itself is 
constant, both the mass and the mass flow rate through the outer parts of the 
active disc decay exponentially on the viscous time scale 



3^'visc(-Rd) 



(8) 



at the outer radius of the active disc. If initially —Md{Rd,t = 0) S> Me, we 
have i?d = RhiM-^) = const, as long as — Md > Mg and the disc is essentially 
emptied during the super-Eddington phase which lasts for a time 



^outb — ''Vise 



-Md(i?d,t = 0) 
Me 



(9) 



This is the behaviour to be expected if the active part of the disc is dominated 
by irradiation, i.e. if i^visc is given by (5). If, on the other hand, viscous heating 
still dominates the vertical temperature structure, z/yisc is given by (4). Using 

(4) in (8) we see that — Md oc M^*'^^, i.e. that the mass flow rate and the mass 
in the active part of the disc do not decay exponentially but rather as 



Mdit) = Md(0) 



1 + 



3^ 
7^ 



-3/7 



(10) 



where 



To 



MM 
-Md(0) 



(11) 



is the time scale on which the disc mass decays initially. If, initially, — Md ^ Me, 
the disc is essentially emptied during the super-Eddington phase which lasts in 
this case for the time 
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where 



touth^ ^{Tm^- To), (12) 



Me 



(13) 



is the time scale on which the disc mass decreases once —M^ reaches Me- 

The duration of the quiescent phase is given by the time it takes to replenish 
the matter which has been lost (accreted/ejected) from the active part of the 
disc during the previous outburst. Because we have assumed that the active disc 
is fed at the rate —M2, irrespective of whether i?disc is larger or smaller than 
Rd, the duration of quiescence is 

t - ^ (14) 
The total length of an outburst cycle is then 

^cycle ~ ^outb ~^ (1^) 

and the duty cycle of the outbursts 

d = ^. (16) 

^cycle 

To decide the question of how much the neutron star can accrete, the quantity 

wc need is the accretion efficiency rj. Because the mass accreted over one cy- 
cle is AMaccr — ^outt-ME, whereas the mass transferred over the same time is 
AMtransf - icycic(-M2), we have 

ioutb-ME ^ Me 



tcycle(-M2) (-M2) 

Writing 



(18) 

Me 



we obtain 



lux 

rj ~ (19) 

X 

if the mass flow rate decays exponentially, i.e. if (5) is used, and 

7 x^-^ - 1 

^-o 03 20 

3x — x^-^ 

if (4), i.e. (10) holds. In both cases, it is immediately seen that 77 is small if 
X » 1, i.e. if during the initial phases of an outburst the mass flow rate is highly 
sup er-Eddington . 
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4. The mass flow rate through the outer disc 

In order to get rj, we need to know x, i.e. —M^{Rd, t = 0). Assuming again that 

Rd = Ri,{Me) < i?disc, -Md{Rd,t = 0) > Me, and using (2), (5), (6) and (8) 
we find 

-Md{Rd,t = 0)= 3.3 X lO-6M0yr-iahac~^^'^ 

\0.6) \6500Kj \"^^^) ^1 

We note that Md{Rd,t = 0) as given in (21), and therefore also rj, still depends 
on the orbital period. The dependence on P enters in a subtle way via the 
filling factor which itself must be a function of M2 and hence via (1) of 
P. Unfortunately, this dependence is not explicitly known and can only be 
determined from full time-dependent calculations of the disc evolution. 

Replacing a^^^^T in (21) by means of (6), noting that 

Md = -M2 • iq, (22) 

we get rid of the poorly known quantities ccc and T at the price of introducing 
fq. This yields 

Md(i?d,t = 0)= 7.6xlO-9MQyr-i(Ce-CR)-Mc^^)"'(6OT) ^23) 
(6imE)-^/^m^'/'m^-74V''''iq(yr). 

Inserting typical values in (23), i.e (Ce — Cr) — l,Q;h — 0.2, rhE = 1.5, mi = 

1.4, m2 = 1 we find x ~ 1.25tq(yr) at P = 20*^ and x ~ 10.6fq(yr) at P = 200'^. 
This shows that if the quiescent time of these discs is long, i.e. > a few decades, 
the accretion efficiency becomes very small. Unfortunately, we do not have 
any direct observational information about tq. However the analogy with the 
properties of the outbursts in black hole soft X-ray transients (BHSXTs) (see 
e.g. King h Ritter 1998) suggests that the quiescent times are likely to be very 
long, i.e. longer than several decades. In the Discussion below wc show that 
the absence of any observed outbursts from the progenitors of pulsars in wide 
circular binaries supports this conclusion. We note in passing that —Ma{Ra) 
must not be too large if this picture is to be self-consistent: the assumption that 
z^visc is dominated by external irradiation, i.e. that we may use (5), requires that 

-^d(-Kh(ME)) < 3fj^g^ ^2^^ 



< 



1.2 X lO-^Moyr-i (g^) (61 mE)2/2mr^ 



If instead of (5) wc use (4), i.e. a viscosity where the disc temperature is self- 
consistently determined by viscous dissipation, we obtain: 

-Md(i?d, i = 0) = 2.0 X lO-^Moyr-i (Ce - UY^'"'' a^'' (^.^ 

(61 mE)-25/28,„-5/^^m|4893pl-3259tq(yr)10/r. ^ > 
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Inserting the same typical values as above in (25) we obtain x ~ 0.7tq(yr)^°/''' 
at P = 20'^ and x ~ 15tq(yr)^°/'^ at P = 200^^. From this it is seen that if only 
tq ^few yr, disc accretion becomes even more inefficient if the viscosity is due 
to viscous dissipation rather than external irradiation. We also note that even 
if (24) is violated initially, the disc's evolution will eventually be dominated by 
external heating because — Md(-Rh) decreases with time. 

Finally, in the case of irradiation-dominated viscosity (5) the duty cycle 
becomes 

d ^ ^visc(Rh(ME)) 

^ 0.095«i7i i^y' (6imE)VV"iqHyr)lnx. ^^^^ 

Although at a first glance d does not seem to depend on P it nevertheless does: 
besides a weak dependence on P via In x there is a significant one from the factor 
tq~^ = —M2/M(i. On the one hand M2 depends explicitly on P (cf (1)), on the 
other hand a more subtle dependence which we have mentioned already earlier 
enters via the filling factor !F in (6). Nevertheless, if typical parameters are used 
in (26), i.e. ~ 0.2, rn-E = 1-5 and mi = 1.4, we have d ~ 0.12 t~^{yT) ■ \nx. 
Therefore, if tq ~ many decades, d <C 1. 

5. Discussion and Conclusions 

Li & Wang (1998) have shown that the accretion efficiency of neutron stars 
which accrete from a giant is very small and, in fact, becomes too small for 
spinning up the neutron star to millisecond spin periods in systems with a final 
orbital period P{ ^ 100'^ if the duty cycle d ^ 10^^. Here we have worked out 
explicit expressions for d and the accretion efficiency 77 in terms of the binary 
parameters and the quiescent time tq of the outburst cycle in the framework of 
an irradiated a-disc model. Our calculations show that both, d and rj are small 
if P > Pcrit and the mass flow rate at the onset of an outburst is highly super- 
Eddington. Whether or not d and rj are really small enough to prevent neutron 
stars in systems with Pf ^ 200*^ being spun up to ms-periods depends crucially 
on the duration of quiescence tq. If tq ^ many decades, as some observations 
of BHSXTs suggest, then both r] and d are very small indeed. Unfortunately, 
determining tq in the framework of our simple model is not possible. For this 
fully time-dependent calculations of the disc instability limit cycle woTild be 
needed. 

However there is another way that we can get a lower limit on d. We know 
that pulsars in wide circular binaries must have descended from SXTs. Yet we do 
not know of a single transient with a period longer than 11.8*^ (GRO J1744-28), 
even though X-ray satellites would probably have detected an outburst from 
such a system anywhere in the Galaxy within the last ~ 30 yr. This must mean 
that outbursts are rather infrequent, i.e. tq is long. We can make this more 
precise by defining the following quantities: 

Let n be the frequency of X-ray outbursts from all systems in the Galaxy, 
A^PSR,obs the currently known number of the pulsar binaries in question, A^psR 
the current number of the pulsar binaries in question in the Galaxy, Tpj-og the 
lifetime of a progenitor system in the phase of nuclear time scale mass transfer. 
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and rpsR the lifetime of the pulsar in the pulsar binaries in question. Here Tprog 
and rpsR are suitable averages over the relevant population. 

Now, the inventory of the binary pulsars in question is incomplete because of 
flux limitation, dispersion and beaming, whereas an X-ray outburst in one of the 
progenitor systems would probably be seen wherever in the Galaxy the outburst 
occurs. Wc correct for this incompletncss by introducing a 'filling factor' / < 1. 
In a stationary situation we then must have 

, _ -^PSR Vog _ -^PSR.obs Vog 

Iq — —J. ■ {■^1 ) 

n rpsR n / trsr 

rpsR must be of the order of the spin-down time scale P/2P ^ 10^ yr, Tprog, 
on the other hand, is of the order of the nuclear time scale toe defined in Ritter 
(1999, Eq. 14) which, for the systems in question, i.e. those which end their 
evolution with a long orbital period Porb ^ 200*^, is of the order 10^ yr. With 
A^PSR,obs = 2 (i.e. PSR B0820+02 and PSR J1803-2712), / ^ 0.2 and from the 
fact that we have not seen a single outburst in the past 30 yr this yields the lower 
limit tq > (60 yr)// ^ 300 yr. We would also like to emphasize that we have 
adopted a very conservative value for /. The real value is likely to be smaller still 
and thus the estimate for t^ even larger because already the canonical beaming 
factor for pulsars is of order 0.2. 

Thus it is very likely that the lack of ms-pulsars in binaries with white 
dwarf companions and Porb ^ 200*^ is entirely a consequence of dwarf nova-like 
disc instabilities in long-period binaries. 

Finally, we note that spinning up a NS in a long-period binary faces yet 
another obstacle: during the long quiescent phases the NS acts as a propeller 
(lUarionov Sz Sunyaev 1975) and thereby is spun down, making the spin-up to 
ms spin periods even more difficult. 
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